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ABBREVIATIONS 


DAMA 

Demand Assignment Multiple Access 

FOV 

Field Of View 

MA 

Multiple-Access 
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New Mexico State University 
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Tracking and Data Relay Satellite 

WSC 
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BACKGROUND 


Currently, the NASA’s Space Network (SN) nominal user services are pre- 
scheduled to allocate service type and required equipment to a user of the network. The 
availability of a resource for a given user is then communicated back to the user prior to 
the requested service time. The current system does not easily allow a real time changes 
to add new users to the service schedule. Therefore a proposed Demand Assignment 
Multiple Access (DAMA) system for NASA’s SN to improve the existing SN usage 
mode is underway at New Mexico State University (NMSU). 

The goals of this project are to improve the current operational modes by 

(a) - reducing the support overhead required to maintain the overall 

scheduling requirements in this area; 

(b) - allowing for more users to access the system especially those who 

traditionally have not considered SN access available to them; and 

(c) - providing additional scheduling flexibility. 

The document “Space Network Demand Assignment Techniques”, prepared by 
Dr. Stephen Horan, Center for Space Telemetering and Telecommunications Systems at 
New Mexico State University, March 16 1995, presents a baseline concept for how this 
DAMA structure might be accomplished. 
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TDRS MA PHASED ARRAY ANTENNA SIMULATIONS 


INTRODUCTION 


One of the baseline concepts discussed in [1] is the S-band Multiple 
Access (SMA) capabilities of the SN. The current configuration for the return MA 
service is to use it in a spot beacon mode as shown in the left half of Figure 1 - Current 
and Proposed MA Coverage for DAMA Request. The user ephemeris is used to predict 
the position of the spot beacon and to track the movement of the user spacecraft. 



Figure 1 - Current and Proposed MA Coverage for DAMA Request 

In the return service signal processing electronics at the White Sands Complex 
(WSC) ground station in New Mexico, weighting factors are added to the return signal 
from each SMA antenna element to provide the correct phase addition to form the spot 
beam. One set of weighting factors and associated signal processing equipment is 
required for each user being serviced by a Tracking and Data Relay Satellite (TDRS) 
within the SN. For the DAMA system, it is proposed to configure the SMA phased array 
to provide a global beacon from each TDRS in this SN as shown in the right half of 
Figure 1. This would require two sets of beam forming hardware at the WSC: one for 
each the east and west TDRS to provide a total of two global beacons with a zone of 
exclusion over the Indian Ocean as is found in the SA services. Two questions arise: Can 
the MA beam formers or antenna elements be configured to provide this global beam and 
if so, what are the performance characteristics of that beam? This report essentially 
answers these two questions by analysing and simulating the TDRS MA Antenna 
therefore producing different beam patterns by varying the effect of the elements on the 
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phased array antenna. It is to be kept in mind that more investigation has to be performed 
to finalise this work as stipulated in the recommendation section. Nonetheless, the basic 
idea is to come out with a beam pattern that has a minimum of loss in the Field Of View 
(FOV) of the TDRS. 

First, some theory will be explained on the TDRS user geometry, the TDRS MA 
antenna and how the simulation program was written. Afterwards the results of the 
simulations of the MA antenna, conclusions and recommendations will be given. 

THEORY 


To be able to simulate the beam pattern of the TDRS Phased- Array Antenna, a 
program was written in Matlab Version 4.2b for Windows. Before showing the program 
and its results, some basic notions which were used to write this program must be 
understood. 

TDRS User Geometry 

Figure 2, shows the geometry or representative TDRS-User Ranges. 


R°fov = 13.5° 



Figure 2. TDRS-User Ranges 


where the earth radius (in degrees) seen by the spacecraft is calculated in the following 
maimer: 


R EARTH — S ^ n 


J 6318KM 
\ 42113 KM - 


= 8.7° « 8.5° 
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Also, the Filed Of View (FOV) of the spacecraft is given by 


R poy — sin 


. JmiiKM' 


V42173£M. 


= 13.9° « 135° 


The approximation are made since the distances are not exact; these 
approximations are only for simulation purposes. Other variables or geometry that will 
be used for the simulation of this array antenna will be the angles: <j> (phi) and 9 (theta). 
These angles will be defined as: 

<j> = the elevation angle from the Nadir ( boresight); and 
0 = the clockwise rotation angle from + x (east). 

. Graphically these angles are represented and referenced as 

POINT OF INTEREST 


FACE PLANE 
OF 

SPACECRAFT 

y = south 



Figure 3 - <j> and 9 Angles Representation 

As an example, White Sands Ground Terminal (WSGT) is at <j) = 8.312° and 0 = 
215.91° (for the east TDRS). With these angles, it will then be possible to determine a 
direction cosine vector for each elements of the array antenna as it will be seen in the next 
section, the TDRS Array Antenna. 

TDRS Array Antenna 

A TDRS is able to support up to 20 separate return user services simultaneously 
through its multiple-access (MA) antenna. This antenna is a phased array of 30 helical 
antenna elements in fixed positions on the face of the antenna module on the TDRS as 
shown on the following figure (Figure 4- TDRS with MA Antenna). 
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Figure 4. TDRS with MA Antenna 

The signal which is received by each of the 30 elements is then transmitted to the 
ground segment via a separate frequency channel. Processing in the ground segment 
extracts the various user return signals by “forming a beam” utilizing the phase 
relationships resulting the 30-element geometry and the TDRS-to-user geometry 
(explained earlier). The beamforming is done entirely within the ground segment. The 
phase relationship between the 30 elements channel must be known accurately in order to 
calculate the phase shifts and form the receive beam. 

The following figure shows the emplacement of the different elements on the MA 
antenna. As shown the coordinate system is right -handed with +x along the flight vector 
(east), +y pointed south, and +z pointed toward the Earth. Actually these coordinates 
correspond to the roll, pitch, and yaw axes respectively on the satellite. 



Figure 5 - Coordinate System of the MA Elements on the Antenna 
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The next table. Table I - Coordinates of the MA Elements, gives the 
nominal x,y,z positions of the receive elements on the TDRS (east) array antenna, with 
respect to the face plane of the spacecraft. All these values are in inches. 


ANTENNA 

X 

y 

V. 

ELEMENTS 

COORDINATES 

COORDINATES 

COORDINATES 


| 27.000 

-38.966 

0.000 


| 13.500 

-38.966 

0.000 


| -1.500 

-37.966 

0.000 


j -13.500 

-38.966 

0.000 

MiMlfii 

-27.000 

-38.966 

0.000 


33.750 

-27.275 

0.000 


20.250 

-27.275 

0.000 


6.750 

-27.275 

0.000 


-6.750 

-27.275 

0.000 


| -20.250 

-27.275 

0.000 

—Milling 

-33.750 

-27.275 

0.000 


| 40.500 

-15.584 

0.000 


27.000 

-15.584 

0.000 


13.500 

-15.584 

0.000 


0.000 

-15.584 

0.000 


-13.500 

-15.584 

0.000 


-27.000 

-15.584 

0.000 


-40.500 

-15.584 

0.000 


13.500 

-2.084 

■ 0.000 

llgtllilglllg 

1 0.000 

-2.084 

0.000 


-13.500 

-2.804 

0.000 

nglMIlMMlg 

-40.500 

11.500 

0.000 


41.413 

14.000 

0.000 


27.913 

14.000 

0.000 


14.413 

14.000 

0.000 


-35.500 

24.000 

0.000 


34.663 

25.691 

0.000 

^Mgilligiiii 

| 21.163 

25.691 

0.000 


| 27.913 

37.382 

0.000 

30 

| 14.413 

37.382 

0.000 


Table I. Coordinates of the MA Elements 


The file ELNOM.DAT in text format contains these values: see Appendix A p. A- 
5/5. These values will be used to determine the direction of the arriving plane wave to 
the antenna. In fact, using the angle defined in the previous section i.e. <j> (phi) and 9 
(theta), Direction Cosine (DC) vectors can be defined. In vector analysis, a Direction 
Cosine is simply the magnitude of the component of the vector along one of the 
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coordinate axes divided by the absolute magnitude of the vector. The nominal DC were 
derived in [2], 


The following DC for x,y and z are similar to those in [2] but the angles were 
defined differently: 


DC(x ) = sin(<j>)cos(0); 
DC(y ) = sin(<j>)sin(0);ura7 
DC(z) = cos(<)>). 


(equation 1) 


Thus the direction or distances from the signal plane to the elements can be obtained by 
multiplying the position of the elements in respect to the spacecraft with the DC defined 
above: 


d = 

"*i Ti *i" 

X 2 T2 Z 2 

X 

'DC(x)' 
DC(y ) 


x i Ti z i" 

x 2 y 2 z 2 

X 

sin(<j))cos( 0 ) 

sin(<j))sin( 0 ) 


_ X 30 T 30 Z 30_ 


DC(z) 


_ X 30 ^30 Z 30 . 


cos(<))) , 


(equation 2) 


o o 

As an example: if the antenna is pointed toward Nadir then <j) = 0 and 0 = 270 then 


DC(x) = sin(0°) cos(270°) = 0; 
DC(y) = sin (0 ) sin(270 ) = 0; and 
DC(z) = cos(0°)= 1. 


Knowing the distance from the signal plane to the elements using the Direction Cosines, 
the relative time of flight can be calculated as: 



(equation 3) 


where c is the speed of light i.e. 2.9986338x10 s m/s. Since the distance (d) is in inches, 
the speed of light (c) will be converted to inches/nanosec given that 1 inch = 0.0254 m 

then c = 2.9986338x10 s m/s x (1 inch/ 0.0254 m) 
c = 1 1 .805644x 1 0 9 inch/sec 


or c = 1 1 .805644 inch/nanosec. 


Also the phase signal can be calculated using 


© = 27t/ 0 T 


rads 


(equation 4) 
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where t is the time of flight calculated earlier and f 0 is the operating frequency (f 0 = 
2.2875 Ghz will be used as the downlink frequency). 


The weights of each elements can then be calculated as (refer to [2]): 


weights(el) = cos(0) - j sin(0) 


(equation 5) 


where “el” stands for element. To simulate the patterns of these elements the Fourier 
Model and Geometric Model explained in detail in [2] will be used. As specified in these 
models, to generate the gain pattern of one element the “sine” function will be utilized: 


G = (sinc(jc)) 2 = j 


(equation 6) 


where for the helical elements on this antenna 



(equation 7) 


and <|) is the elevation angle in degrees as defined earlier. 


Graphically we have: 



Nominal Pattern 
of one element 


Figure 6 - Antenna Element Pattern 

As an example, it is stipulated in [3] that the maximum gain of one element is 16 dB 
(boresight) and the nominal gain should be down to 13.3 dB at 13.5°. 
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Using equations 6 and 7 above: 


for (j) = 0 (Boresight) 


then x = — 7i = 0 
32 


G 0 = (sinc(O)) 2 


sin(0)V =1 

( 0 ) 


(using Hospital’s Rule) 


for $ = 13.5° 


13i° 

then x = ^=-7t = 1.3254 
32 


13 . 5 “ 


= ( sin ^ 1 - 3254 )) 2 = 


sin(1.3254)V 
(1.3254) , 


= 0.535678 


Finally 

[G 0 ] - [G 13 5 ] = 101og(l) - 101og(0.535678) 

= 0 - (-2.71) 

= 2.71 dB 

which meets the specification (given in [3]) of 

16 dB- 13.3 dB = 2.7dB. 

The 16 dB represents the maximum gain of the main lobe of one element and 
therefore 13.5° is approximately the Half Power Bandwidth (HPBW) since there is a loss 
of « 3 dB. 

Also these results were compared with theory using [4]. If 

C = circumference of the helix; 

S = spacing between turns = Ctan(a); 
a = pitch angle = tan' 1 (S/C); 

L = length of one turn = - Jc 2 + S 2 ; 

N = number of turns; 

A = axial length = NxS; and 
d = diameter of the helix conductor. 
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The following figure shows these different parameters. 




Axis 

of 

helix 


Figure 7. Helical Element Parameters (from reference [4]) 
If a frequency (downlink) of 2.2875 Ghz is used then 


X = — = 


c = f x A. = speed of light 
3x 10 8 m / s 
f 22875 xlO 9 Hz 


= 0.131w 


using some of the practical results from [4]: C = 0.92A. and a = 1 3° 
then C -0.92x0.131 = 0.12 m 


and 

S = Ctan(a) = 0.12xtan(13°) = 0.028 m 

Also since these elements have N = 21 turns and from [4], the Half Power 
Bandwidth (HPBW) is given by: 



(equation 8) 


Then replacing the numerical values in this last expression: 


HPBW = — 

V0.13y V v 0.1317 


HPBW = 26.8° 
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which corresponds approximately to the 13.5° (HPBW/2 = 13.4°) half power bandwidth 
calculated earlier. Also to calculate the gain, let’s first find the directivity using the 
following formula from [4] again: 



(equation 9) 


thus 


D = 15- 



21-0.028 

0.131 


D = 565 


As stipulated in [4], “this directivity expression is a ratio (not decibel), and 
approximately equals the gain since the axial mode helix is nearly lossless”. Then 


G = 10 log (D) 


(equation 10) 


G = 101og(56.5) = \152dB 


which is similar to the specification gain that was mentionned earlier i.e. 16 dB. 
References [5] anf [6] were also used for these calculations. 

SIMULATION OF THE MA ANTENNA 

Using all the theory mentionned in the previous sections and using part of 
reference [3] routines, a program was written in Matlab Version 4.2b for Windows to 
simulate the beam pattern of the MA antenna. This program was called PATTERN.M 
and the listing is in Appendix A page A-2/5 to A-5/5. The following parts of the program 
PATTERN.M must be noted. First this program needs the file ELNOM.DAT which 
contains the positions of the elements of the phased-array antenna with respect to the face 
plane of the spacecraft. Appendix A p. A-5/5 shows the content of this file which 
contains the same numbers presented in Table I in this report. Also the amplitude in dB 
given as the output in the simulations is in reference to the maximum gain of one element 
i.e. 16 dB as discussed earlier in the report. Thus if the maximum amplitude of a pattern 
is -1 dB, this means that it is 1 dB below the reference of 16 dB therefore the real gain is 
15 dB. Table II- Summary of the Simulations (p.20/22) shows the values of HPBW and 
Gain for the different simulations. Finally to simulate for various angle of <|> (phi) and 0 
(theta), these variables must be changed: see Change Variables Section in PATTERN.M. 

A few simulation were performed and they will be briefly explained in the 
following paragraphs. 
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First Simulation: 30 Elements at Nadir ((j) = 0°. 6=270°1 

The first simulation produced the beam pattern of the 30 elements combined 
together when the antenna is pointing at Nadir. To be able to simulate these elements 
some modifications had to be done in the program in Appendix A. First since the 
simulation was going to be performed at the Nadir, the “phi” and “theta” variables had 
to be given values of 0° and 270° respectively as shown in the Change Variables Section 
in PATTERN.M. Also the section EFFECT OF 30 ELEMENTS page A-4/5 was created 
and the variable “iel” in the for loop was identified as iel=l :30. It is to be noted that the 
other sections EFFECT OF ELEMENT # 20, EFFECT OF ELEMENTS # (14-15-16) and 
EFFECT OF ELEMENT # (19-20-21) were put in comments (“%” is the comment 
symbol in Matlab). Appendix B, p.B-2/24 to B-4/24 shows the results of the simulation 
of these elements (at Nadir). The first plot Appendix B p.B-2/22 shows that the gain 
varies significantly depending were we are situated in the FOV. This kind of pattern is 
not what is desired since the gain is not really constant over the FOV. The next plot (p.B- 
3/24) shows a transversal 2-D plot of this simulation to be able to observe the main lobe. 
The maximum gain of this main lobe is approximately 30 dB. The plot p.B-4/24 shows a 
contour plot of this simulation: top view. 

Second Simulation: 30 Elements at WSGT (<|) = 8.312° and 6 = 215.91°) 

This second simulation gives similar results to the first one in respect to the fact 
that the gain is not constant. To be able to simulate these elements, some modifications 
were also done in the program. First, since the simulation was going to be performed at 
different angles than for Nadir, the “phi” and “theta” variables had to be given values of 
8.312° and 215.91° respectively in the Change Variables Section in PATTERN.M. Also 
the section EFFECT OF 30 ELEMENTS page A-4/5 was kept identical as the last 
simulation. It is to be noted that the other sections EFFECT OF ELEMENTS were put 
in comments. Three plots are shown in page B-5/24 and p.B-7/24. 


Third Simulation: Element #20 at Nadir (<j) = 0° and 9 = 270°1 

The third simulation was performed by using only one element and see its effect 
over the FOV. The element chosen was element number 20 (central element in the 
phased-array antenna). First, since the simulation was going to be performed at the 
Nadir, the “phi” and “theta” variables had to be given values of 0° and 270° 
respectively as shown in the Change Variables Section in PATTERN.M. Also the section 
EFFECT OF ELEMENT # 20 page A-4/5 was created and the variable “iel” in the for 
loop was identified as iel=20:20. It is to be noted that the other sections EFFECT OF 
ELEMENTS were put in comments ( “%”). Appendix B, p.B-8/24 to B-l 1/24 show the 
results of the simulation of this element (at Nadir). The first plot shows a 3 Dimensional 
representation of the element over a range of -50° to 50° in the x and y axis. The main 
lobe has a gain of 0 dB relative to the reference of 16 dB as stipulated earlier. Also it was 
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mentionned earlier that the FOV of the TDRS is approximately 13.5° therefore a plot 
focusing only in that region is shown in the next page p.B-9/24. As seen the gain is 
almost flat in that range as desired. The next plot (p.B- 10/24) shows a transversal (2-D) 
cut of this last plot : Gain in dB (z axis) and FOV in degrees (x axis). The loss at 13.5° is 
less than 3 dB (approximately 2.7 dB as calculated earlier in the report). These results are 
symmetric as shown in the next plot (p.B-1 1/24) where z axis is shown versus the y axis. 

Fourth Simulation: Element #20 at WSGT & = 8.312° and 6 = 215.91°! 

The fourth simulation was performed by using element number 20 again. This 
time the simulation was performed at <|) = 8.312° and 0 = 215. 91° which is the 
emplacement of WSGT, therefore the “phi” and “theta” variables had to be given values 
of 8.312° and 215.91° respectively in the Change Variables Section in PATTERN.M. 

All the rest of the program was kept the same as the last simulation. 

The first plot p.B- 12/24 shows a 3 Dimensional representation of the element over 
a range of -50° to 50° in the x and y axis: as it can be seen this plot is very similar to the 
previous simulation at Nadir. The next plot shows a zoom in the -13.5° and 13.5° range 
(p.B- 13/24). As seen the gain is almost flat in that range as desired (very similar to plot 
p.B-9/24). The next two plots show transversal (2-D) cuts of this last plot : Gain in dB(z 
axis) and FOV in degrees (x axis) and FOV in degrees (y axis). Again , these plots are 
very similar to the last simulation i.e. Element # 20 at Nadir 

Fifth Simulation: Elements #19-20-21 at Nadir = 0° and 0 = 270°) 

The fifth simulation was performed by using the 3 central elements of the phased- 
array antenna and see the effect of combining the beam patterns over the FOV. The 
elements that were chosen are elements number 19, 20 and 21. First since the simulation 
was going to be performed at the Nadir, the “phi” and “theta” variables had to be given 
values of 0° and 270° respectively as shown in the Change Variables Section in 
PATTERN.M. Also, the section EFFECT OF ELEMENT #(19-20-21) page A-4/5 was 
created and the variable “iel” in the for loop was identified as iel= 19:21. The other 
sections EFFECT OF ELEMENTS were put in comments ( “%”). Appendix B, p.B- 
16/24 and B- 17/24 show the results of the simulation of these elements combined 
together (at Nadir). The first plot shows a 3 dimensional (3-D) representation of the 
element over a range of -13.5° to 13.5° (desired FOV) in the x and y axis. The main lobe 
has a gain of approximately 9.5 dB relative to the reference of 16 dB as stipulated earlier. 
Thus the actual gain is approximately 25.5 dB (9.5 dB +16 dB). As seen the gain is not 
flat in that range (-13.5° to 13.5°) therefore the effect of adding other elements tends to 
produce sidelobes. The next plot shows a transversal (2-D) cut of this last plot : Gain in 
dB(z axis) and FOV in degrees (x axis). 
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Sixth Simulation: Elements #19-20-21 at WSGT (4> = 8.312° and ft = 215.91°^ 


The sixth simulation also was performed using the 3 central elements of the 
phased-array antenna as the fifth simulation but this time the simulation was performed at 
WSGT. Since the simulation was going to be performed at the WSGT, the “phi” and 
“theta” variables had to be given values of 8.312° and 215.91° respectively in the 
Change Variables Section in PATTERN.M. Appendix B, p.B-18/24 and B-19/24 show 
the results of the simulation of these elements combined together (at WSGT). The first 
plot shows a 3 Dimensional representation of the element over a range of -13.5° to 13.5° 
(desired FOV) in the x and y axis. This plot is quite different from the last simulation at 
Nadir. It seems that the main lobe has been shifted. Again, the gain is not flat in that 
range which shows the effect of adding other elements tends to produce sidelobes. The 
next plot p.B- 19/24 shows a transversal (2-D) cut of this last plot : Gain in dB(z axis) and 
FOV in degrees (x axis). The maximum gain of the main lobe is approximately 9.5 dB in 
reference to 16 dB again. 

Seventh Simulation: Elements #(14-15-161 and (19-20-21) at Nadir t(j) = 0° and 6 = 270°1 

This simulation was performed by using the 6 central elements of the phased- 
array antenna and see the effect of combining the beam patterns over the FOV. The 
elements that were chosen are elements number 14, 15, 16, 19, 20 and 21. Again, since 
the simulation was going to be performed at the Nadir, the “phi” and “theta” variables 
had to be given values of 0° and 270° respectively as shown in the Change Variables 
Section in PATTERN.M. Also the sections EFFECT OF ELEMENT #(14-15-16) and 
EFFECT OF ELEMENT #(19-20-21) page A-4/5 were used and the variable “id” in the 
for loop were identified as iel=14:16 and iel= 19:21 respectively. The other sections 
EFFECT OF ELEMENTS were put in comments ( “%”). Appendix B, p.B-20/24 and B- 
21/24 show the results of the simulation of these elements combined together (at Nadir). 
The first plot shows a 3 dimensional representation of the elements over a range of -13.5° 
to 13.5° (desired FOV) in the x and y axis. The main lobe has approximately a gain of 

15.5 dB relative to the reference of 16 dB as stipulated earlier (total gain: approximately 

3 1 .5 dB). Again, as seen on the plot the gain is not flat in that range therefore the effect 
of adding other elements tends to produce sidelobes as dicussed earlier. But it is to be 
noted, that the symmetry is better than with 3 elements only. The next plot, p.B-21/24, 
shows a transversal (2-D) cut of this last plot : Gain in dB(z axis) and FOV in degrees (x 
axis). 

Eight Simulation: Elements #(14-15-16) and (19-20-21) at WSGT 
t<|) = 8.132° and 9 = 215.91° 

The last simulation also was performed by using the 6 central elements of the 
phased-array antenna as the seventh simulation but this time the simulation was 
performed at WSGT. Since the simulation was going to be performed at the WSGT, the 
“phi” and “theta” variables had to be given values of 8.312° and 215.91° respectively as 
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shown in the Change Variables Section in PATTERN.M. Appendix B, p.B-22/24 to B- 
24/24 show the results of the simulation of these elements combined together (at WSGT). 
The first plot shows a 3 Dimensional representation of the element over a range of -13.5° 
to 13.5° (desired FOV) in the x and y axis. This plot is quite different from the last 
simulation at Nadir. It seems that the main lobe has been shifted. Again, the gain is not 
flat in that range therefore the effect of adding other elements tends to produce sidelobes. 
The next two plots p.B-23/24 and B-24/24 show transversal (2-D) cuts of this last plot : 
Gain in dB(z axis) versus FOV in degrees (x axis) and FOV in degrees (y axis) 
respectively. 

The following table, Table II. Summary of the Simulations, gives a comparison of 
the main lobes of the different simulation. 


MAIN LOBE • 1 

SIMULATION 

RELATIVE 

MAXIMUM 

GAIN 

HPBW/2 

GAIN at 
HPBW 
(* 3dB less) 

COMMENTS 

FIRST SIMULATION 
(30 elements at Nadir) 
see p.B-2/24 to p.B-4/24 

29.5 dB 

* 1° 

26.5 dB 

gain is not flat in 
the desired FOV i.e 
-13.5° to 13.5° 

SECOND SIMULATION 
(30 elements at WSGT) 
see p.B-5/24 to p.B-7/24 

29.5 dB 

* 1° 

26.5 dB 

gain is not flat in 
the desired FOV i.e 
-13.5° to 13.5° 

THIRD SIMULATION 
(Element #20 at Nadir) 
see p.B-8/24 to p.B-1 1/24 

0 dB 

« 13.5° 

-2.7 dB 

gain is flat in FOV 
but the relative 
max gain is lower 

FOURTH SIMULATION 
(Element #20 at WSGT) 
see p.B-1 2/24 to p.B-1 5/24 

*1.5 dB 

* 13.5° 

*- 1.5 dB 

gain is flat in FOV 
but the relative 
max gain is lower 

FIFTH SIMULATION 
(Elements #19-20-21 at Nadir) 
see p.B-1 6/24 and p.B-1 7/24 

9.5 dB 

depends on what axis 
x axis : « 1° 
y axis : * 13.5° 

* 6.5 dB 

gain not flat in the 
desired FOV 

SIXTH SIMULATION 
(Elements #19-20-21 at WSGT) 
see p.B-1 8/24 and p.B-1 9/24 

9.5 dB 

depends on what axis 
x axis : * 1 ° 
y axis : * 13.5° 

* 6.5 dB 

main lobe is shifted 
compared to the 
fifth simulation 

SEVENTH SIMULATION 
(Elements #(14-15-16) & #(19- 
20-21) at Nadir) 
see p.B-20/24 and p.B-21/24 

15.5 dB 

« 3.5° 

12.5 dB 

gain not flat in the 
desired FOV 

EIGHTH SIMULATION 
(Elements #(14-15-16) & #(19- 
20-21) at WSGT) 
see p.B-22/24 to p.B-24/24 

15.5 dB 

« 3.5° 

12.5 dB 

main lobe is shifted 
compared to the 
seventh simulation 


Table II. Summary of the Simulations 

Note : Table II is referenced to 16 dB (gain relative to a single element). Thus the true 
value of the gains in the table are found by adding 16 dB to any gain. 
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CONCLUSION 


It seems that throughout all the simulations performed and shown in this report, 
the best result for a whole coverage over the FOV of the TDRS is obtained by using only 
one element (in this case element # 20 was simulated). Nonetheless, it is to be noted that 
the simulation of elements #(14-15-16) & (19-20-21) gave good results even if it does not 
cover the whole FOV. Certainly by performing simulations choosing different 
combinations of elements, results similar but not identical to what was found with one 
element would be acquired. The results will be similar and not identical since by adding 
more elements, sidelobes are created in the FOV as noted in the simulations. 


RECOMMENDATION 

• The TDRS MA Antenna should be looked into more closely to see if only one 
element has enough power to cover the whole FOV. A link budget should therefore 
be performed; 

• More simulations should be done to be able to determine the best pattern covering the 
whole FOV with a minimum of attenuation; and 

• Actual measurements with the antenna should be accomplished to compare the real 
measurements of gain and FOV to the simulations in this report; 
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APPENDIX A 


PROGRAM 


Listing of program PATTERN.M p.A-2/5 

Position of the MA elements on the TDRS antenna: 

file ELNOM.DAT p.A-5/5 


A- 1/5 



o / 0 *** ********************************************************** * 
0 ^***** *************************** ****************************** 
0/o ******** PROGRAM : pattem.m ****** 

%******** SOFTWARE USED: Matlab Version 4.2b for Windows ****** 

%******** DATE : 10 April 1995 ****** 

%******** DESCRIPTION : This program will simulate the TDRS **** 

%* ****** * phased-array antenna elements * * * * 

0 ^***** ******************************** ****************** ******* 
o/ Q ************************************************************** 


clear; % clearing all variables in the Matlab environment 

load elnom.dat; % file containing the NOMinal ELement positions ... 
% ...on the MA TDRS antenna 

% CHANGE VARIABLES SECTION 

phi= 0; % elevation angle from Nadir 

theta = 270; % clockwise rotation angle from +x (east) 

% INITIALIZE VARIABLES 

be = 1 ; % index position of x,y and z for the antenna... 

iy = 2; % ... elements. 

iz = 3; 


phi_min =-13.5; % minimum value for the elevation angle phi 

phijmax =13.5; % maximum value for the elevation angle phi 

phi inc = 0.5; % increment for the angle phi 

floor_gain = -30; 

dataline = zeros( 1 ,4 096); % initialize the gain vector 

% USE THE VALUES OF PHI AND THETA TO PRODUCE THE BEAM PATTERN 


xsinc = phi/32; 

gnom = 10 * loglO((sinc(xsinc)) A 2); 
rphi = pi*phi/180; 
rtheta = pi*theta/180; 
dc(ix) = cos(rtheta)*sin(rphi); 
dc(iy) = sin(rtheta)*sin(rphi); 
dc(iz) = cos(rphi); 


% equation 7 page 13/22 in report 
% equation 6 page 13/22 in report 
% convert phi in degrees to radians 
% convert theta in degrees to radians 
% Direction Cosine equations... 

% ... see equation 1 page 12/22 ... 

% ... in report 
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% CALCULATE THE WEIGHTS OF EACH ELEMENTS 
% (based on the conjugate of the desired signal) 

for iel= 1:30 

disnom = dc(ix)*elnom(iel,ix)+dc(iy)*elnom(iel,iy)+dc(iz)*einom(ie!,iz); 

timenom = disnom/1 1 .805644; % 1 1 .805644 is speed of light in inches/nanosec 

% ... timenom is the relative time of flight... 

% ... to the element: see report page 12/22 ... 

% ... equation 3. 

phsnom = 2*pi*timenom*2.2875; % phsnom is the phase signal where 2.2875 GHz... 

% ... represents the downlink frequency: ... 

% ... see report page 12/22 equation 4. 

weights(iel)= cos(pihsnom)- j*sin(phsnom); 

% weight of each element: equation 5 ... 

% ... report page 13/22. 

end; %for iel 

% CREATE THE BEAM PATTERN 
column =0; 

for xaxis = phi_min:phi_inc:phi_max % x axis pattern 

column = column + 1; 
dataline = xaxis; 
lineinc =0; 
yaxis = phimin; 

while ((yaxis >= phi min) & (yaxis <= phi_max)) % y axis pattern 
loop = 0; 

phi = sqrt(xaxis A 2 + yaxis A 2); 

while ((phi > phi max) & (yaxis <= phi_max)) 
line_inc = Iine_inc + 1 ; 
dataline(Iine_inc)= floor_gain; 
yaxis = yaxis + phiinc; 
phi = sqrt(xaxis A 2 + yaxis A 2); 
loop = loop + 1 ; % we were in the loop 

end % while phi & yaxis 

if(phi = 0) 
theta = 0; 
else 

theta = atan2(yaxis,xaxis); % the result of atan2(y,x)... 

% ... is in the interval ... 

% ...[pi -pi]. 

end %if phi=0 

xsinc = phi/32; 
rphi = pi*phi/1 80; 

dc(ix) = cos( theta)* sin(rphi); % calculate the DC for the ... 

dc(iy) = sin(theta)*sin(rphi); %... values of phi and theta... 
dc(iz) = cos(rphi); %... to plot the beam pattern, 

signom = 0; 
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% EFFECT OF ELEMENT #20 
% for iel = 20:20 

% disnom = dc(ix)*elnom(iel,ix)+ dc(iy)*elnom(iel,iy)+ dc(iz)*elnom(iel,iz); 

% timenom = disnom/1 1 .805644; 

% phsnom = 2*pi*timenom*2.2875; 

% signom = signotn + (cos(phsnom)+ j*sin(phsnom))* weights(iel); 

% end; %for iei = 20 

% EFFECT OF ELEMENTS # (14-15-16) 

%for iel = 14:16 

% disnom = dc(ix)*elnom(iel,ix)+dc(iy)*elnom(iei,iy)+dc(iz)*elnom(iel,iz); 

% timenom = disnom/1 1 .805644; 

% phsnom = 2*pi*timenom *2.2875; 

% signom = signom + (cos(phsnom)+ j*sin(phsnom))* weights(iel); 

%end; %for iel 14 to 16 

% EFFECT OF ELEMENTS # (19-20-21) 

%for iel = 19:21 

% disnom = dc(ix)*elnom(iel,ix)+dc(iy)*elnom(iel,iy)+dc(iz)*elnom(iel,iz); 

% timenom = disnom/1 1 .805644; 

% phsnom = 2*pi*timenom*2.2875; 

% signom = signom + (cos(phsnom)+ j*sin(phsnom))* weights(iel); 

%end; %for iel 19 to 21 

% EFFECT OF 30 ELEMENTS 
for iel =1:30 

disnom = dc(ix)*elnom(iel,ix)+dc(iy)*eInom(iel,iy)+dc(iz)*elnom(iel,iz); 
timenom = disnom/1 1 .805644; 
phsnom = 2 *pi*timenom *2.2875; 

signom = signom + (cos(phsnom)+ j*sin(phsnom))* weights(iel); 
end; %for iel 1 to 30 

num_el = 1 ; % number of elements 

gain = 10*logl0(((abs(signom)*sinc(xsinc)) A 2)/(num_el) A 2)-gnom; % total gain 

if (gain < floor_gain) 
gain = floor_gain; 
end; %if gain 

if loop = 0 

lineinc = iine_inc + 1; 
dataline(lineinc) = gain; 
yaxis = yaxis + phiinc; 
end; %if loop = 0 

end; %while yaxis 

xx(l:(line_inc), column) = (ones(l,Hne_inc)*xaxis)'; % setting x axis values to plot 
yy(l:(line_inc), column) = (phi_min:phi_inc:phi_max)'; % setting y axis values to plot 
zz( I :(line_inc), column) = (dataline)’; % setting z axis values to plot 


end; % for xaxis 

mesh(xx,yy,zz); % plot the Beam Pattern in 3-D 

% plot(elnom(:,l),elnom(:,2),'*'); % plot the MA Antenna Phased-Array Element Placement 
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FILE: ELNOM.DAT 


% Description : This file contains the coordinates x,y,z of the MA elements on the TDRS 
Phased- Array Antenna with respect to the center of the spacecraft. 


27.0000 

-38.966 0 

13.5000 

-38.966 0 

-1.500 • 

-37.966 0 

-13.500 

-38.966 0 

-27.000 

-38.966 0 

33.750 

-27.275 0 

20.250 

-27.275 0 

6.750 - 

27.275 0 

-6.750 ■ 

-27.275 0 

-20250 

-27.275 0 

-33.750 

-27.275 0 

40.500 

-15.584 0 

27.000 

-15.584 0 

13.500 

-15.584 0 

0.000 - 

15.584 0 

-13.500 

-15.584 0 

-27.000 

-15.5840 

-40.500 

-15.584 0 

13.500 

-2.084 0 

0.000 

-2.084 0 

-13.5000 -2.084 0 

-40.500 

11.500 0 

41.413 

14.000 0 

27.913 

14.000 0 

14.413 

14.000 0 

-35.500 

24.000 0 

34.663 

25.691 0 

21.163 

25.691 0 

27.913 

37.382 0 

14.413 

37.382 0 
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APPENDIX B 


MATLAB BEAM PATTERN SIMULATIONS 

- First Simulation: 

30 Elements at Nadir (<J> = 0°, 0=270°) p.B-2/24 

- Second Simulation: 

30 Elements at WSGT (<j> = 8.312° and 0 = 215.91°) p.B-5/24 

- Third Simulation: 

Element #20 at Nadir (<j) = 0° and 0 = 270°) p.B-8/24 

- Fourth Simulation: 

Element #20 at WSGT (<j> = 8.312° and 0 = 215.91°) p.B-12/24 

- Fifth Simulation: 

Elements #19-20-21 at Nadir (<|> = 0° and 0 = 270°) p.B-16/24 

- Sixth Simulation: 

Elements #19-20-21 at WSGT (<j> = 8.312° and 0 = 215.91°) p.B-18/24 

- Seventh Simulation: 

Elements #(14-15-16) and (19-20-21) 

at Nadir (<f> = 0° and 0 « 270°) p.B-20/24 

- Eight Simulation: 

Elements #(14-15-16) and (19-20-21) 

at WSGT (<j> = 8.312° and 0 = 215.91°) p.B-22/24 


B-l/24 



Gain in dB (z axis) 


FTRST SIMULATION 


30 Elements at Nadir (phi=0 degrees and theta=270 degrees) 
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Gain in dB (z axis) 
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Gain in dB (z axis) 



30 Elements at WSGT (phi=8.312 degrees and theta=215.91 degrees) 



B-5/24 


SECOND SIMULATION 


30 Elements at WSGT (phi=8.312 degrees and theta=215.91 degrees) 
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THIRD SIMULATION 


Element #20 at Nadir (phi=0 degrees and theta=270 degrees) 



FOV in degrees (y axis) 


-50 


-50 


FOV in degrees (x axis) 



THIRD SIMULATION 


Element #20 at Nadir (phi=0 degrees and theta=270 degrees) 



FOV in degrees (x axis) 
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THIRD SIMULATION 


Element #20 at Nadir (phi=Q degrees and theta=27Q degrees) 



FOV in degrees (y axis) 
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FOURTH SIMULATION 


Element #20 at WSMR (phi=8.312 degrees and theta=215.91 degrees) 
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FOURTH STMIJLATTON 


Element #20 at WSMR (phi=8.312 degrees and theta=215.91 degrees) 



FOV in degrees (y axis) 


-20 -15 


-10 


FOV in degrees (x axis) 
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FOURTH SIMULATION 


Element #20 at WSMR (phi=8.312 degrees and theta=215.91 degrees) 
5r 



FOV in degrees (x axis) 
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FOURTH SIMULATION 


Element #20 at WSMR (phi=8.312 degrees and theta=215.91 degrees) 
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Gain in dB (z axis) 



Elements #19-20-21 at Nadir (phi= 0 degrees and theta = 270 degrees) 
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FTFTH SIMULATION 


Elements #19-20-21 at Nadir (phi= 0 degrees and theta = 270 degrees) 



FOV in degrees (x axis) 
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Gain in dB (z axis) 


Elements #19-20-21 at WSGT (phi= 8.312 degrees and theta = 215.91 degrees) 
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Gain in dB (z axis) 



B- 19/24 


Gain in dB (z axis) 


Elements #(14-15-16) & (19-20-21) at Nadir (phi=0 degrees and theta=270 degrees) 
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SEVENTH SIMULATION 


Elements #(14-15-16) & (19-20-21) at Nadir (phi=0 degrees and theta=270 degrees) 



FOV in degrees (x axis) 
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Gain in dB (z axis) 



Elements #(14-15-16) & (19-20-21) at WSGT (phi=8.132 degrees and theta=215.91 degrees) 
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Gain in dB (zaxis) 
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EIGHT SIMULATION 


Elements #(14-15-16) & (19-20-21) at WSGT (phi=8.132 degrees and theta=215.91 degrees) 
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